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Abstract. Present cosmological observations yield an upper bound on the neutrino mass which is sig- 
nificantly stronger than laboratory bounds. However, the exact value of the cosmological bound is model 
dependent and therefore less robust. Here, I review the current status of cosmological neutrino mass bounds 
and also discuss implications for sterile neutrinos and LSND in particular. 
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1 Introduction 



The absolute value of neutrino masses are very difficult 
to measure experimentally. On the other hand, mass dif- 
ferences between neutrino mass eigenstates, (mi, m2, m.3), 
can be measured in neutrino oscillation experiments. 

The combination of all currently available data sug- 
gests two important mass differences in the neutrino mass 
hierarchy. The solar mass difference of <5m 12 — 7 x 10 -5 
eV 2 and the atmospheric mass difference Sm^ — 2.6 x 
1CT 3 eV 2 [THUS] (see also the contribution by C. Giunti 
to the present volume). 

In the simplest case where neutrino masses are hierar- 
chical these results suggest that mi <~ 0, m 2 ~ <5m so i ar , 
and m.3 ~ i5m a t m ospheric- If the hierarchy is inverted 0flSJ 
IqTT7II8TI9"] one instead finds m 3 ~ 0, m 2 ~ <5m atmospheric , 
and mi ~ (Jmatmospheric- However, it is also possible that 
neutrino masses are degenerate |10II11II12II13II14II15II16II171 
IT%||19II20|. mi ~ m 2 ~ m 3 > <5m atm os P heric, hi which case 



oscillation experiments are not useful for determining the 
absolute mass scale. 

Experiments which rely on kinematical effects of the 
neutrino mass offer the strongest probe of this overall mass 
scale. Tritium decay measurements have been able to put 
an upper limit on the electron neutrino mass of 2.2-2.3 eV 
(95% conf.) [2] (see also the contribution by C. Kraus 
in the present volume). However, cosmology at present 
yields an even stronger limit which is also based on the 
kinematics of neutrino mass. 

Neutrinos decouple at a temperature of 1-2 MeV in the 
early universe, shortly before electron-positron annihila- 
tion. Therefore their temperature is lower than the pho- 
ton temperature by a factor (4/11) 1 / 3 . This again means 
that the total neutrino number density is related to the 



photon number density by 



11 



(1) 



Massive neutrinos with masses m ^> To ~ 2.4 x 1CP 4 
eV are non-relativistic at present and therefore contribute 
to the cosmological matter density 22,23,24 



n v h 2 = 



92.5 eV 



(2) 
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calculated for a present day photon temperature T = 
2.728K. Here, J2 m v = m i + m 2 + m 3 . However, because 
they are so light these neutrinos free stream on a scale 

of roughly k ~ 0.03m cV ^m /2 h Mpc -1 |2oTI^6ll27| . Below 
this scale neutrino perturbations are completely erased 
and therefore the matter power spectrum is suppressed, 
roughly by AP/P ~ -%n v /Q m g7|. 

This power spectrum suppression allows for a deter- 
mination of the neutrino mass from measurements of the 
matter power spectrum on large scales. This matter spec- 
trum is related to the galaxy correlation spectrum mea- 
sured in large scale structure (LSS) surveys via the bias 
parameter, b 2 = P g (k)/P m (k). Such analyses have been 
performed several times before |28II29| . most recently us- 
ing data from the 2dF galaxy survey |3()j . 

However, using large scale structure data alone does 
not allow for a precise determination of neutrino masses, 
because the power spectrum suppression can also be caused 
by changes in other parameters, such as the matter density 
or the Hubble parameter. 

Therefore it is necessary to add information on other 
parameters from the cosmic microwave background (CMB). 
This has been done in the past 30 31 32 , using ealier 
CMB data. More recently the precise data from WMAP 
[53] has been used for this purpose [4TJIl34ll35| to derive a 
limit of 0.7-1.0 eV for the sum of neutrino masses. 
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2 Cosmological data and likelihood analysis 

The extraction of cosmological parameters from cosmo- 
logical data is a difficult process since for both CMB and 
LSS the power spectra depend on a plethora of differ- 
ent parameters. Furthermore, since the CMB and matter 
power spectra depend on many different parameters one 
might worry that an analysis which is too restricted in 
parameter space could give spuriously strong limits on a 
given parameter. 

The most recent cosmological data is in excellent agree- 
ment with a flat /1CDM model, the only non-standard fea- 
ture being the apparently very high optical depth to reion- 
ization. Therefore the natural benchmark against which 
non-standard neutrino physics can be tested is a model 
with the following free parameters: i7 m , the matter den- 
sity, the curvature parameter, J?b, the baryon density, Ho, 
the Hubble parameter, n s , the scalar spectral index of 
the primordial fluctuation spectrum, r, the optical depth 
to reionization, Q, the normalization of the CMB power 
spectrum, b, the bias parameter, and finally the two pa- 
rameters related to neutrino physics, Q v h 2 and N v . The 
analysis can be restricted to geometrically flat models, i.e. 
fl = fi m + fl/y = 1. For the purpose of actual power 
spectrum calculations, the CMBFAST package can 
be used. 



2.1 LSS data 

At present, by far the largest survey available is the 2dF- 
GRS 37 of which about 147,000 galaxies have so far been 
analyzed. Tegmark, Hamilton and Xu [3H1 have calculated 
a power spectrum, P(k), from this data, which we use 
in the present work. The 2dFGRS data extends to very 
small scales where there are large effects of non-linearity. 
Since we only calculate linear power spectra, we use (in 
accordance with standard procedure) only data on scales 
larger than k = 0.2ft, Mpc -1 , where effects of non-linearity 
should be minimal 39 . Making this cut reduces the num- 
ber of power spectrum data points to 18. 



2.2 CMB data 

The CMB temperature fluctuations are conveniently de- 
scribed in terms of the spherical harmonics power spec- 
trum Ci = (|a im | 2 ), where 4f(6,<f>) = Y^i m a lm Y lm{0,4>)- 
Since Thomson scattering polarizes light there are also 
power spectra coming from the polarization. The polar- 
ization can be divided into a curl-free (E) and a curl 
(B) component, yielding four independent power spectra: 
Ct,i, Ce,Ii Cb,i and the temperature i?-polarization cross- 
correlation Cte,i- 

The WMAP experiment have reported data only on 
C T ,i and Ctej, as described in Ref. j3^1l4Tl4Tll4^1l4^j 

We have performed the likelihood analysis using the 
prescription given by the WMAP collaboration which in- 
cludes the correlation between different Cj's j53T Hrlll4*Tl 



Table 1. 95% C.L. upper limits on for the three 

different cases: 1) WMAP+Wang+2dFGRS+HST+SN-Ia, 2) 
WMAP+Wang+2dFGRS 3) WMAP+2dFGRS. 



Case 


(95% C.L.) 


1 


1.01 eV 


2 


1.20 eV 


3 


2.12 eV 



42 43]. Foreground contamination has already been sub- 
tracted from their published data. 

In parts of the data analysis we also add other CMB 
data from the compilation by Wang et al. |44j which in- 
cludes data at high I. Altogether this data set has 28 data 
points. 



3 Neutrino mass bounds 

The analysis presented here was originally published in 
Ref. [H3]> and more details can be found there. 

We have calculated % 2 as a function of neutrino mass 
while marginalizing over all other cosmological parame- 
ters. This has been done using the data sets described 
above. In the first case we have calculated the constraint 
using the WMAP C T ,i combined with the 2dFGRS data, 
and in the second case we have added the polarization 
measurement from WMAP. Finally we have added the ad- 
ditional constraint from the HST key project and the Su- 
pernova Cosmology Project. It should also be noted that 
when constraining the neutrino mass it has in all cases 
been assumed that N v is equal to the standard model 
value of 3.04. Later we relax this condition in order to 
study the LSND bound. 

The result is shown in Fig. 1. As can be seen from 
the figure the 95% confidence upper limit on the sum of 
neutrino masses is Y2, m v < 1.01 eV (95% conf.) using the 
case with priors. This value is completely consistent with 
what is found in Ref. j35j where simple Gaussian priors 
from WMAP were added to the 2dFGRS data analysis. 
For the three cases studied the upper limits on "^2 m v can 
be found in Table 1. 

In the middle panel of Fig. 1 we show the best fit value 
of Hq for a given fi u h 2 . It is clear that an increasing value 
of m v can be compensated by a decrease in Hq. Even 
though the data yields a strong constraint on f2 m h 2 there 
is no independent constraint on fl m in itself. Therefore, 
an decreasing Hq leads to an increasing fl m . This can be 
seen in the bottom panel of Fig. 1. 

When the HST prior on Hq is relaxed a higher value 
of m v is allowed, in the case with only WMAP and 
2dFGRS data the upper bound is Q v h 2 < 0.023 (95% 
conf.), corresponding to a neutrino mass of 0.71 eV for 
each of the three neutrinos. 
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This effect was also found by Elgar0y and Lahav in 
their analysis of the effects of priors on the determination 
of ^m„. 

However, as can also be seen from the figure, the ad- 
dition of high-Z CMB data from the Want et al. compila- 
tion also shrinks the allowed range of m v significantly. 
The reason is that there is a significant overlap of the 
scales probed by high-Z CMB experiments and the 2dF- 
GRS survey. Therefore, even though we use bias as a free 
fitting parameter, it is strongly constrained by the fact 
that the CMB and 2dFGRS data essentially cover much 
of the same range in /c-space. 

It should be noted that Elgar0y and Lahav [23 find 
that bias does not play any role in determining the bound 
on ^2 m v . At first this seems to contradict the discussion 
here, and also what was found from a Fisher matrix analy- 
sis in Ref. [31] . The reason is that in Ref. [3S] , redshift dis- 
tortions are included in the 2dFGRS data analysis. Given 
a constraint on the amplitude of fluctuations from CMB 
data, and a constraint on fl m h 2 , this effectively constrains 
the bias parameter. Therefore adding a further constraint 
on bias in their analysis does not change the results. 

Neutrinoless double beta decay - Recently it was claimed 
that the Heidelberg-Moscow experiment yields positive ev- 
idence for neutrinoless double beta decay. Such experi- 
ments probe the 'effective electron neutrino mass m ee = 
| Ugjirivj |. Given the uncertainties in the involved nu- 
clear matrix elements the Heidelberg-Moscow result leads 
to a mass of m ee = 0.3—1.4 eV. If this is true then the mass 
eigenstates are necessarily degenerate, and Y m v — 3m ee . 
Taking the WMAP result of ]T m v < 0.70 eV at face value 
seems to be inconsistent with the Heidelberg-Moscow re- 
sult 0H]- However, already if Ly-a forest data and a con- 
straint on the bias parameter is not used in the analysis 
the upper bound of Y, rn f < 1.01 eV is still consistent. 
For this reason it is probably premature to rule out the 
claimed evidence for neutrinoless double beta decay. 

Evidence for a non-zero neutrino mass - In a recent 
paper |4fi| it was noted that there is a preference for a non- 
zero neutrino mass if a measurement of the bias parameter 
from X-ray clusters is added to the CMB and large scale 
structure data. This result arises because the X-ray data 
prefers a low value of as (bias), which is incompatible with 
the WMAP and 2dF result at the 2a level. While this is 
an interesting finding it is clear that the X-ray data is 
subject to a serious problem with systematic uncertainties, 
such as the calibration of the mass-temperature relation. 
Therefore the result more likely points to a problem with 
the interpretation of the X-ray data than to evidence of a 
non-zero neutrino mass. 



4 Sterile neutrinos 

In Ref. [j2] it was shown that there is a degeneracy be- 
tween the neutrino mass (X] m ^) an d the relativistic en- 
ergy density, parameterized in terms of the effective num- 
ber of neutrino species, N v . 

As can be seen from Fig. 2, the best fit actually is ac- 
tually shifted to higher Y m v when N v increases, and the 
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Fig. 1. The top panel shows x 2 as a function of ^Jn„ for 
different choices of priors. The dotted line is for WMAP + 
2dFGRS data alone, the dashed line is with the additional 
Wang et al. data. The full line is for additional HST and SNI- 
a priors as discussed in the text. The horizontal lines show 
A\ 2 = 1 and 4 respectively. The middle panel shows the best 
fit values of Ho for a given m v ■ The horizontal lines show 
the HST key project la limit of H = 72 ± 8 kmsMpc" 1 . 
Finally, the lower panel shows best fit values of C2 m . In this 
case the horizontal line corresponds to the SNI-a la upper 
limit of Sl m < 0.42. 



conclusion is that a model with high neutrino mass and 
additional relativistic energy density can provide accept- 
able fits to the data. As a function of N v the upper bound 
on Yj m v ( a t 95% confidence) can be seen in Table 2. 

This has significant implications for attempts to con- 
strain the LSND experiment using the present cosmologi- 
cal data. Pierce and Murayama conclude from the present 
MAP limit that the LSND result is excluded 021 ( see a l so 
Ref. 0Z|). 

However, for several reasons this conclusion does not 
follow trivially from the present data. In general the three 
mass differences implied by Solar, atmospheric and the 
LSND neutrino measurements can be arranged into either 
2+2 or 3+1 schemes. Recent analyses 0S] of experimental 
data have shown that the 2+2 models are ruled out. The 
3+1 scheme with a single massive state, J714, which makes 
up the LSND mass gap, is still marginally allowed in a 
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Table 2. 95% C.L. upper limits on ^2,m v for different values 
of N v . 



effective N„ 


J2m v (95% C.L.) 


3 


f .01 eV 


4 


1.38 eV 


5 


2.12 eV 



few small windows in the {Am 2 , sin 2 29) plane. These gaps 
are at {Am 2 , sin 2 20) ~ (0.8 eV 2 , 2 x 1CT 3 ), (1.8 eV 2 , 8 x 
l(T 4 ),(6eV 2 ,l-5 x 1CT 3 ) and (10eV 2 ,1.5 x lO" 3 ). These 
four windows corresponds to masses of 0.9, 1.4, 2.5 and 3.2 
eV respectively. From the Solar and atmospheric neutrino 
results the three light mass eigenstates contribute only 
about 0.1 eV of mass if they are hierarchical. This means 
that the sum of all mass eigenstate is close to m^. 

The limit for N v = 4 which corresponds roughly to the 
LSND scenario is ^m u < 1.4 eV, which still leaves the 
lowest of the remaining windows. The second window at 
m ~ 1.8 eV is disfavoured by the data, but not at very 
high significance. 




0.0050.010.0150.020.0250.03 



Fig. 2. Ax 2 as a function of ^2 m v for various different values 
of N v . The full line is for N„ = 3, the dotted for N u = 4, and 
the dashed for N„ = 5. A\ 2 is calculated relative to the best 
fit N v = 3 model. 



5 Discussion 

We have calculated improved constraints on neutrino masses 
and the cosmological relativistic energy density, using the 



new WMAP data together with data from the 2dFGRS 
galaxy survey. 

Using CMB and LSS data together with a prior from 
the HST key project on H yielded an upper bound of 
J2 m u < 1-01 eV (95% conf.). While this excludes most of 
the parameter range suggested by the claimed evidence for 
neutrinoless double beta decay in the Heidelberg-Moscow 
experiment, it seems premature to rule out this claim 
based on cosmological observations. 

Another issue where the cosmological upper bound on 
neutrino masses is very important is for the prospects of 
directly measuring neutrino masses in tritium endpoint 
measurements. The successor to the Mainz experiment, 
KATRIN, is designed to measure an electron neutrino 
mass of roughly 0.2 eV, or in terms of the sum of neu- 
trino mass eigenstates, X) m ^ — 0-75 eV (see contribution 
by Guido Drexlin to the present volume). The WMAP 
result of m v — 0.7 eV (95% conf.) already seems to ex- 
clude a positive measurement of mass in KATRIN. How- 
ever, this very tight limit depends on priors, as well as 
Ly-a forest data, and the more conservative present limit 
of J2 m u — 1-01 e V (95% conf.) does not exclude that 
KATRIN will detect a neutrino mass. 

Finally, we also found that the neutrino mass bound 
depends on the total number of light neutrino species. In 
scenarios with sterile neutrinos this is an important factor. 
For instance in 3+1 models the mass bound increases from 
1.0 eV to 1.4 eV, meaning that the LSND result is not 
ruled out by cosmological observations yet. 
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